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The process pp —> 7+jet + ^ is studied using 1.0 fb^^ of data collected by the DO detector at the 
Fermilab Tevatron pp collider at a center-of-mass energy =1.96 TeV. Photons are reconstructed in 
the central rapidity region < 1.0 with transverse momenta in the range 30<Py <400 GeV while 
jets are reconstructed in either the central ly""'! < 0.8 or forward 1.5 < ly-""'! < 2.5 rapidity intervals 
with pJy' > 15 GeV. The differential cross section d^a/dp^dy^dy^^^ is measured as a function of 
in four regions, differing by the relative orientations of the photon and the jet in rapidity. Ratios 
between the differential cross sections in each region are also presented. Next-to-leading order QCD 
predictions using different parameterizations of parton distribution functions and theoretical scale 
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choices are compared to the data. The predictions do not simultaneously describe the measured 
normalization and pip dependence of the cross section in the four measured regions. 

PACS numbers: 13.85.Qk, 12.38.Qk 



The production of a photon with associated jets in 
the final state is a powerful probe of the dynamics of 
hard QCD interactions [l|, 0, S, 0, H H ■ Different angu- 
lar configurations between the photon and the jets can 
be used to extend inclusive photon production measure- 
ments 0, 0, H [IBl and simultaneously test the underlying 
dynamics of QCD hard-scattering subprocesscs in differ- 
ent regions of parton momentum fraction x and large 
hard-scattering scales . 

In this Letter, we present an analysis of photon plus 
jets production in pp collisions at a center-of-mass en- 
ergy ^/s =1.96 TeV in which the most-energetic (leading) 
jhoton is produced centrally with a rapidity \y'^\ < 1.0 



xi < 0.024, 0.110 < X2 < 



0.300 for y'^-yj'^' 



> and 



11| . The cross section as a function of photon trans- 
verse momentum p^ is measured differentially for four 
separate angular configurations of the highest px (lead- 
ing) jet and the leading photon rapidities. The leading 
jet is required to be in either the central {\y^°^\ < 0.8) 
or forward (1.5 < \y^°*'\ < 2.5) rapidity intervals, with 
> 15 GeV, and the four angular configurations stud- 

> and with < 0, 



ict 

Ft 

icd are: central jets with y''-y^^^ 
and forward jets with y^ ■y^'^'^ 



> and with y'^-yj^' < 0. 
The total x and region covered by the measurement is 
0-007 < X < 0.8 and 900 < = {p'^.f < 1.6 x 10^ GeV^, 
extending the kinematic reach of previous photon plus jet 

HIilHIil- Ratios between 



12, 13, 14 



measurements 

the differential cross sections in the four studied angular 
configurations are also presented. The measurements are 
compared to the corresponding theoretical predictions. 

Isolated final-state photons produced in pp — > 7-l-jet-l- 
X events are expected to mainly originate "directly" from 
QCD Compton-likc qg qj scattering or qq 57 an- 
nihilation subprocesscs. In Fig. [T] the expected contri- 
bution, estimated using pythia [I3 Monte Carlo (MC) 
event generator with the CTEQ6.5M parton distribution 
function (PDF) set [2^, of the Compton-like partonic 
scattering process to the total associated production of 
a photon and a jet is shown for each of the four mea- 
sured rapidity intervals. The parton distribution func- 
tions entering into the theoretical predictions have sub- 
stantial uncertainties, particularly for the gluon contri- 
butions at small x, large x and large 2^, 21 1. The 
measurement intervals probe different regions of parton 
momentum-fraction space of the two initial interacting 
partons, xi^2- For example at = 40 GeV, in events 
with a central leading jet, the y^ -y^^^ > region covers 
adjacent xi and X2 intervals (0.016 ^ < 0.040 and 
0.040 <X2< 0.100), while for events with y'^ ■y^''^ < 0, 
the xi and X2 intervals are similar (0.029 ^ a;i < 0.074, 
0.027 ^ 2^2 ^ 0.065). In events with a forward leading 
jet, intervals of small and large x are covered (0.009 < 



0-097 < XI < 0.264, 0.022 < X2 < 0.059 for y'^ -y^"^ < 0). 
Here xi^2 are defined using the leading order approxima- 



tion Xi^2 = (Pt/\/^)('^ 



±y-' 



') [ii,aii,sj]- 



The data presented here correspond to an integrated 
luminosity of 1.01 ± 0.06 fb"! ^ collected using the DO 
detector at the Fermilab Tevatron pp collider operating 
at a center-of-mass energy ^/s =1.96 TeV. A detailed 
description of the DO detector can be found in [i^l and 
only an overview of the detector components relevant to 
this analysis is given here. 

Photon candidates are formed from clusters of 
calorimeter cells in the central part of the liquid-argon 
and uranium calorimeter. The central calorimeter cov- 
ers the pseudorapidity range \r]\ < 1.1 and two end 
calorimeters cover 1.5 < \ri\ < 4.2 The electro- 

magnetic (EM) section of the central calorimeter con- 
tains four longitudinal layers of 2, 2, 7, and 10 radiation 
lengths, and is finely-segmented transversely into cells of 
size Ar]X A(j> ^ 0.1 x 0.1 (0.05 x 0.05 in the third EM 
layer), providing good angular resolution for photons and 
electrons. The position and width of the Z boson mass 
peak, reconstructed from Z — > e~^e~ events, are used 
to determine the EM calorimeter calibration factors and 
the EM energy resolution [2^ . The central section of the 
calorimeter surrounds a central preshower detector, with 
three concentric cylindrical layers of scintillator strips, 
and a tracking system consisting of silicon microstrip and 
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FIG. 1: The fraction of events, estimated using the pythia 
event generator, produced via the qg — > q'y subprocess relative 
to the total associated production of a direct photon and a 
jet for each of the four measured configurations of the leading 
jet and leading photon rapidities. 
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FIG. 2: Normalized distribution of the ANN output Onn for 
data, 7 -f jet signal MC, and dijet background MC events for 
44 < < 50 GeV after application of the main selection 
criteria. 
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FIG. 3: Normalized distribution of the ANN output Onn for 
electrons from boson decays in data and MC events. 



scintillating fiber trackers located within a 2 T solcnoidal 
magnetic field. 

The DO tracking system is used to select events which 
contain a primary collision vertex, reconstructed with at 
least three tracks, within 50 cm of the center of the de- 
tector along the beam axis. The efHciency of the vertex 
requirement varies as a function of instantaneous lumi- 
nosity from 92% to 96%. 

Photon candidates with rapidity \y^\ < 1.0 are selected 
from clusters of calorimeter cells within a cone of ra- 
dius Tl = y/ (A?7)^ -I- (A(/))2 = 0.4 defined around a seed 
tower [i^l • The final cluster energy is then re-calculated 
from the inner cone with TZ = 0.2. The data are se- 
lected using a combination of triggers based on photon 
EM shower profiles in the calorimeter and EM cluster px 
thresholds. The total trigger efficiency is (96-97)% for 
photon candidates with « 32 GeV and greater than 
99% for > 40 GeV. The selected clusters are required 
to have greater than 96% of their total energy contained 
in the EM calorimeter layers. Isolated clusters are se- 
lected by requiring that the energy Eem{T^ = 0.2), calcu- 
lated within the inner cone of radius TZ = 0.2, fulfills the 
condition [Etote^iin = 0.4) - E^Min - 0.2)]/£;EM(7^ = 
0.2) < 0.07, where £;total(7^ = 0.4) is the summed EM 
and hadronic energy within a cone of radius TZ = 0.4. 
The candidate EM cluster is required not to be spatially 
matched to a reconstructed track. This is accomplished 
by computing a function evaluating the consistency, 
within uncertainties, between the reconstructed 77 and 
(f> positions of the cluster and the closest track in the 
finely-segmented third layer of the EM calorimeter. The 
corresponding probability is required to be < 0.1%. 
Background contributions to the direct photon sample 
from cosmic rays and from isolated electrons, originat- 
ing from the leptonic decays of W bosons, important 



at high [2^ , are suppressed by requiring the missing 
transverse energy ^t, calculated as a vector sum of the 
transverse energies of all calorimeter cells, to satisfy the 
condition < 12.54-0.36 pi^ GeV. The longitudinal seg- 
mentation of the EM calorimeter and central prcshowcr 
detector allow us to estimate the photon candidate direc- 
tion and vertex coordinate along the beam axis (" photon 
vertex pointing" ) . This vertex is required to lie within 
10 cm of the event primary vertex reconstructed from 
charged particles. 

Photons arising from decays of tt*^ and ij mesons are 
already largely suppressed by the requirements above, 
and especially by photon isolation, since these mesons 
are produced mainly within jets during fragmentation 
and are surrounded by other particles. To better select 
photons and estimate the residual background, an ar- 
tificial neural network (ANN) is constructed using the 
JETNET package The following three variables are 

used in the ANN: the number of cells in the first EM 
layer belonging to the cluster, the fraction of the cluster 
energy deposited in the first EM layer, and the scalar 
sum of charged particle transverse momenta in the hol- 
low cone 0.05 < TZ < 0.4 around the photon cluster di- 
rection. The resulting ANN output, Onn, after applying 
all data selection criteria, is shown, normalized to unit 
area, in Fig. [2] for 44 < pj, < 50 GeV. The output is 
compared to photon signal events and dijet background 
events simulated using pythia. The signal events may 
contain photons originating from the parton-to-photon 
fragmentation process. For this reason, the background 
events, produced with QCD processes in pythia, were 
preselected to exclude the bremsstrahlung photons pro- 
duced from partons. Signal and background MC events 
were processed through a GEANT-based simulation of 
the DO detector and the same reconstruction code as used 
for the data. The ANN is tested using electrons from 
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FIG. 4: The purity of the selected pp 
and jet rapidities. The resuhs of the 1 - 
uncertainties (dotted hues), and the total uncertainties (dashed lines). 



7 + jst + X sample as a function of for each measured configuration of photon 
exp(a + functional fits are shown by the solid lines, together with the systematic 



Z boson decays and the resulting normalized data and 
MC distributions are compared in Fig. [3l Photon candi- 
dates are selected by the requirement Onn > 0.7 which 
has good background rejection and a signal efficiency in 
the range (93-97)%. The signal selection efficiency de- 
creases by about 4% with increasing from 30 GeV to 
300 GeV due to the Onn > 0.7 requirement. The to- 
tal photon-|-jet selection efficiency after applying all the 
selection criteria, including the ANN and the require- 
ments, is (63-77)% as a function of with an overall 
systematic uncertainty of (4.7-5.2)%. Main sources of 
inefficiency are the isolation, anti-track matching, ANN, 
and the photon vertex pointing cuts. 

Events containing at least one hadronic jet arc selected. 
Jets are reconstructed using the DO Run II jet-finding al- 
gorithm with a cone of radius 0.7 [2^ . and are required to 
satisfy quality criteria which suppress background from 
leptons, photons, and detector noise effects. Jet energies 
are corrected to the particle level. The leading jet should 
havepi^* > 15 GeV and lyj''*] < 0.8 or 1.5 < lyj''*] < 2.5. The 
leading photon candidate and the leading jet are also re- 
quired to be separated in rj—cf) space by A7?,(7, jet) > 0.7. 
The leading jet total selection efficiency varies from 94% 
to almost 100% and takes into account any migrations 
between leading and second jet from the particle to the 
reconstruction level. The total systematic uncertainty 
on this efficiency is 5.7% at pj, ~ 30 GeV, decreasing 
to about 2% at p'^ > 200 GcV. The measurement is not 



very sensitive to jet energy scale corrections since it is 
performed in bins of (with > 30 GeV) and only 
information on the the jet angular direction is used. 

In total, about 1.4 million candidate events are selected 
after application of all selection criteria. A correction for 
the "7-1-jet" event purity V is then applied to account for 
the remaining background in the region Onn > 0.7. The 
distribution of the ANN output for the simulated pho- 
ton signal and dijet background samples are fitted to the 
data for each bin using a maximum likelihood fit (soj 
to obtain the fractions of signal and background com- 
ponents in the data without constraining the fractions of 
signal and background samples in the fit to be in the [0, 1] 
range. The data and fitted sum of the weighted signal 
and background MC distributions of Onn are found to 
be compatible with /ndf values in the range 0.2-1.3 
3l| . The resulting purities are shown in Fig. 2] for each 
measurement region. The dependence of the purity 
is fitted in each region using a two parameter function 
V = 1 — exp(a + bpTp). The result of the fits together 
with their statistical errors are shown in Fig. 21 The 
systematic uncertainties on the fit are estimated using 
alternative fitting functions and varying the number of 
bins in the fitting of the ANN output distribution. An 
additional systematic uncertainty due to the fragmen- 
tation model implemented in PYTHIA is also taken into 
account. It is found to be 5% at ~ 30 GeV, 2% at 
~ 50 GeV, and 1% at pJ. > 70 GeV 0. 
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To study whether bremsstrahlung photons have dif- 
ferent selection efficiencies from direct photons, we ex- 
tracted them from dijet events simulated with PYTHIA. 
We found that they do not produce a noticeable change 
of the selection efficiencies, acceptance and shape of the 
photon ANN output that have been obtained with direct 
photons. 

The differential cross section d"^(T/dpJdj/'''dyj°* for the 
process 7 -f jet + X is obtained from the num- 

ber of data events in each interval, after applying correc- 
tions for background, efficiency, and acceptance effects, 
divided by the integrated luminosity and the widths of 
the interval in the photon transverse momentum, photon 
rapidity, and jet rapidity. The data are also corrected for 
Pt bin-migration effects which result from the finite en- 
ergy resolution of the EM calorimeter using an analytical 
Ansatz method [s^ and the measured EM energy resolu- 
tion determined from the Z boson peak. The correction 
factors range from (1-5)% with about a 1% uncertainty. 

The total {Sa°^^) and main sources of experimental 
systematic uncertainty are shown for the \y^'^*'\ < 0.8, 
yi . yjct ^ Q i-ggiQj^ Fig. [51 Similar uncertainties are 
found for the other measured regions. The largest uncer- 
tainties are assigned to the purity estimation [(10-4)%], 
photon and jet selections [(7.7-5.2)%], photon energy 
scale [(4.2-6.0)%], and the integrated luminosity (6.1%). 
The uncertainty ranges above are quoted with uncer- 
tainty at lowp^- first and at highp^ second. The system- 
atic uncertainty on the photon selection is due mainly to 
the anti-track match cut (3%), a correction due to ob- 
served data/MC difference in the efficiency of the main 
photon selection criteria found from Z —tee events [(1-5- 
2)%], the photon vertex pointing requirement (2%), the 
ANN cut (2%), and the uncertainty on the parameterized 
photon selection efficiency (<1%). The total experimen- 
tal systematic uncertainty for each data point is obtained 
by adding all the individual contributions in quadrature. 

The result for each region is presented as a function 
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FIG. 5: The total and main sources of systematic uncertainty 
for the cross section measured in the < 0.8, y''-y'°^ > 
rapidity region. 
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FIG. 6: The measured differential pp —> 'y + jet -I- X cross 
section as a function of for the four measured rapidity 
intervals. For presentation purposes, the cross section results 



\y 



< 0.8) jets with y'' -y'"^ > and for forward 



< 2.5) jets with y'' ■ y'"^ > and y'' ■ y'"^ < 



for central 
(1.5 < \y>^'\ 

are scaled by factors of 5, 0.1 and 0.3, respectively. The data 
are compared to the theoretical NLO QCD predictions u sing 
the JETPHOX package [H with the CTEQ6.5M PDF set ^ 
and renormalization, factorization and fragmentation scales 
f^R = fiF ^ fj-f =plf{y*)- 



TABLE I: Differential cross sections d^a/dp'^dy'^dy'^^ and un- 
certainties for the <0.8, y'^-y''^^ > rapidity interval. 



p ' 


bin 


w> 


Cross section 


I^Crstat 






(GeV) 


(GoV) 


(pb/GoV) 


(%) 


(%) 


(%) 


30 - 


34 


31.9 


3.08x10^ 


0.2 


14.2 


14.2 


34 - 


39 


36.3 


1.74x10^ 


0.3 


13.1 


13.1 


39 - 


44 


41.3 


9.76x10° 


0.4 


12.4 


12.4 


44 - 


50 


46.8 


5.60x10° 


0.5 


11.9 


11.9 


50 - 


60 


54.6 


2.76x10° 


0.6 


11.5 


11.5 


60 - 


70 


64.6 


1.24x10° 


0.9 


11.0 


11.0 


70 - 


80 


74.7 


6.25 xlO'"- 


1.2 


10.8 


10.9 


80 - 


90 


84.7 


3.32 xlO"^ 


1.7 


10.6 


10.7 


90 - 


110 


99.0 


1.51 xlO"^ 


1.8 


10.6 


10.7 


110 - 


130 


119.1 


5.79 xlO"^ 


2.9 


10.5 


10.9 


130 - 


150 


139.2 


2.56 xlO"^ 


4.3 


10.7 


11.5 


150 - 


170 


159.3 


1.17 xlO"^ 


6.5 


10.9 


12.7 


170 - 


200 


183.6 


5.80 xlO"^ 


7.6 


11.0 


13.3 


200 - 


230 


213.8 


2.33 xlO"^ 


11.8 


11.0 


16.1 


230 - 


300 


259.5 


7.25 xlO"" 


13.8 


10.7 


17.5 


300 - 


400 


340.5 


7.96 xlO~^ 


35.3 


10.9 


36.9 



of in Fig. [H and Tables [TlUVl The data points are 
plotted at the value (p!^) for which a value of the smooth 
function describing the cross section equals the average 
cross section in the bin [s^] ■ The data cover six orders of 
magnitude in the cross section for events with < 0.8, 
falling more rapidly over four orders of magnitude for 



events with 1.5 < < 2.5 
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FIG. 7: The ratios of the measured triple-differential cross section, in each measured interval, to the NLO QCD prediction using 
JETPHOX Q with the CTEQ6.5M PDF set and all three scales /^a.f,/ ~ P^fiv*)- The solid vertical line on the points shows 
the statistical and pr-dependent systematic uncertainties added in quadrature, while the internal line shows the statistical 
uncertainty. The two dotted lines represent the effect of varying the theoretical scales by a factor of two. The shaded region is 
the CTEQ6.5M PDF uncertainty. The dashed and dash-dotted hues show ratios of the jetphox predictions with MRST 2004, 
Alekhin, and ZEUS 2005 to CTEQ6.5M PDF sets. Systematic uncertainties have large (> 80%) pip bin-to-bin correlations. 
There is a common 7.8% normalization uncertainty that is not shown on the data points. 



The data are compared to next-to-leading order (NLO) 
QCD predictions obtained using jetphox [H, [sll, with 
CTEQ6.5M PDF [2^ and BFG fragmentation func- 
tions of partons to photons 



37| 



The rcnormalization, 
factorization, and fragmentation scales (/ii?, /iF, and 
/^/) are set equal to p^fiu*), where f{y*)^ {[1 
exp 



(-2|y*|)]/2}i/2 and y* = O.biy"' - i/j<=*) [381]. The 
theoretical predictions include selection criteria on the 
photon and jet similar to those applied in the experi- 
mental analysis. In particular, an isolation requirement 
on the photon of [Etotaii'R- = 0.4) - E"']/E'y < 0.07 is 
made, where Etotai{T^ = 0.4) is the total energy around 
the photon in a cone of radius TZ = 0.4, and E'^ is the 
photon energy. This requirement suppresses the relative 
contribution from photons produced in the fragmentation 
process, and leads to a more consistent comparison with 
the experimental result. Corrections for the underlying 



event and parton-to-hadron fragmentation contributions, 
estimated using pythia, are found to be negligibly small 
and are not included. To make a more detailed compar- 
ison, the ratio of the measured cross section to the NLO 
QCD prediction is taken in each interval and the results 
arc shown in Fig. [T] The inner error bars reflect the 
statistical uncertainty only, and the outer error bars are 
the total statistical and pT-dependent systematic uncer- 
tainties summed in quadrature. Most of these systematic 
uncertainties, associated with the paramctrizations of the 
photon and jet selection efhciencies, purity (including 
the uncertainty from the PYTHIA fragmentation model), 
photon pt correction, and calorimeter energy scale, have 
large (> 80%) bin-to-bin correlations in p^. Systematic 
p!^-independcnt uncertainties from the luminosity mea- 
surement, photon selection efficiency caused by the anti- 
track matching. ANN and photon vertex pointing, ac- 
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TABLE II: 



Differential 

\y 



cross sections d'^cr/dpydj/^dj/^°* and 
\y'°^\ <0.8, y'^-y'"^ < rapidity interval. 



Pt 


bin 


(Pt) 


Cross section 


(^CTstat 






(GcV) 


(GoV) 


(pb/GoV) 


(%) 


(%) 


(%) 


30 - 


34 


31.9 


2.51x10^ 


0.3 


15.7 


15.7 


34 - 


39 


36.3 


1.42x10^ 


0.3 


13.9 


13.9 


39 - 


44 


41.3 


7.90x10" 


0.4 


12.6 


12.6 


44 - 


50 


46.8 


4.48x10" 


0.5 


11.9 


11.9 


50 - 


60 


54.6 


2.20x10° 


0.6 


11.5 


11.5 


60 - 


70 


64.6 


9.99x10"^ 


0.9 


11.1 


11.1 


70 - 


80 


74.7 


4.98x10"^ 


1.3 


10.9 


11.0 


80 - 


90 


84.7 


2.67x10"^ 


1.8 


10.7 


10.9 


90 - 


110 


99.0 


1.26x10"^ 


1.9 


10.7 


10.9 


110 - 


130 


119.1 


4.74x10"^ 


3.1 


10.6 


11.1 


130 - 


150 


139.2 


2.07x10"^ 


4.7 


10.9 


11.9 


150 - 


170 


159.3 


1.08x10"^ 


6.6 


11.2 


13.0 


170 - 


200 


183.6 


5.23x10"^ 


7.7 


11.7 


14.0 


200 - 


230 


213.8 


1.90x10"^ 


13.0 


11.6 


17.4 


230 - 


300 


259.5 


5.93x10"'* 


15.0 


11.2 


18.7 


300 - 


400 


340.5 


5.32x10"^ 


46.1 


12.9 


47.8 



TABLE III: Differential cross sections d^cr/dp^^dy^dj/j''* and 



uncertainties for the 
interval. 



1.5 < < 2.5 



p ' 


bin 




Cross section 








(GcV) 


(GoV) 


(pb/GoV) 


(%) 


(%) 


(%) 


30 - 


34 


31.9 


1.67x10' 


0.3 


14.7 


14.7 


34 - 


39 


36.3 


8.74x10° 


0.4 


13.5 


13.5 


39 - 


44 


41.3 


4.53x10° 


0.5 


12.8 


12.8 


44 - 


50 


46.8 


2.36x10° 


0.7 


12.4 


12.4 


50 - 


60 


54.5 


1.02x10° 


0.8 


11.8 


11.8 


60 - 


70 


64.6 


3.96x10"* 


1.4 


11.2 


11.3 


70 - 


80 


74.6 


1.71x10"* 


2.1 


10.8 


11.0 


80 - 


90 


84.7 


7.76x10"^ 


3.2 


10.8 


11.3 


90 - 


110 


98.8 


3.05x10"^ 


3.6 


10.7 


11.3 


110 - 


130 


118.9 


8.27x10"^ 


6.9 


11.0 


13.0 


130 - 


150 


139.0 


2.85x10"^ 


11.8 


11.5 


16.5 


150 - 


200 


169.4 


3.15x10"'* 


23.0 


12.1 


26.0 



ceptance (1.5%), and unfolding (1%) lead to a total 7.8% 
overall normalization uncertainty and are not sho'wn in 

Fig.m 

The prediction using the CTEQ6.5M PDF and BGF 
fragmentation sets does not describe the shape of the 
cross section over the whole measured range. In par- 
ticular, the prediction is above the data for events with 
lyj''*! < 0.8 in the region > 100 GeV and below the 
data for jets produced in the 1.5 < lyj*^*! < 2.5, y'^-y^^'^ > 
rapidity region for < 50 GeV. Most of the data points 

TABLE IV: Differential cross sections d^tj/dp^dy^dj/*''* and 



uncertainties for the 
interval. 



1.5 < < 2.5, y^-y^"' < rapidity 



P ' 


bin 


{pi-} 


Cross section 




ifayat 




(GcV) 


(GoV) 


(pb/GoV) 


(%) 


(%) 


(%) 


30 - 


34 


31.9 


8.08x10° 


0.4 


15.6 


15.6 


34 - 


39 


36.3 


4.36x10° 


0.4 


14.2 


14.2 


39 - 


44 


41.3 


2.23x10° 


0.6 


13.0 


13.0 


44 - 


50 


46.8 


1.16x10° 


0.8 


12.3 


12.3 


50 - 


60 


54.5 


5.28x10"* 


1.0 


11.7 


11.7 


60 - 


70 


64.6 


2.08x10"* 


1.7 


11.3 


11.4 


70 - 


80 


74.6 


9.18x10"^ 


2.6 


11.2 


11.5 


80 - 


90 


84.7 


4.61x10"^ 


3.7 


11.3 


11.9 


90 - 


110 


98.8 


1.64x10"^ 


4.5 


11.2 


12.1 


110 - 


130 


118.9 


5.31x10"^ 


8.2 


11.1 


13.8 


130 - 


150 


139.0 


1.79x10"^ 


14.1 


11.2 


18.0 


150 - 


200 


169.4 


3.04x10"'* 


23.0 


11.3 


25.6 



in these and rapidity regions are (1-1.5) Satot outside 
of the CTEQ6.5M PDF set uncertainty range which is 
shown by the shaded region in the figure and calculated 
according to the prescription in . Note that the data- 
to-theory ratios have a shape similar to those observed 
in the inclusivephoton cross sections measured by the 
UA2 01, CDF [il and DO [l^ cohaborations. 

The dotted Hues in Fig.[7]show the effect of setting the 
renormalization, factorization, and fragmentation scales 
to 0.5p^/(j/*) (upper dotted line) and 2p7pf{y*) (lower 
dotted line). The effect on the normalization is (9-11)%, 
except for jets in the 1.5 < < 2.5, y'^ -y^"* < ra- 
pidity range where it is (18-20)%. The scale variation 
is not able to simultaneously accommodate the mea- 
sured differential cross sections in all of the measured 
regions. The ratios of the NLO QCD prediction with the 
MRST 2004 [ii, Alekhin 0, and ZEUS 2005 [13 PDF 
sets to the prediction obtained using the CTEQ6.5M 
PDF set are also presented in the figure. The shapes 
of the predictions are very similar, especially for forward 
jet production, with the different PDF sets. 

The ratios of the predicted cross sections with the de- 
fault scales [/Xfl = /xf = A*/ = Pt/(2/*)] to those with all 
the scales set equal to p^ are presented for each of the 
four kinematic regions as a function of in Fig. [51 For 
each measured region, the new prediction is smaller than 
the default case across the entire p^ range, most notably 
in the forward jet rapidity intervals where this choice of 
scale leads to a poorer level of agreement between data 
and theory. 

Uncertainties related to the photon production due to 
the fragmentation mechanism are also studied separately 
using the JETPHOX package. The ratio of the pp ^ -f + 
jet + X cross section for the direct photon contribution 
to the sum of direct and fragmentation contributions is 
shown, for the chosen photon isolation criteria, in each 
of the four measured regions in Fig. (S] For all regions, 
the fragmentation contribution decreases with increasing 
pI di, 52, m and is largest for the 1.5 < \y^°^\ < 2.5, 
yf . yjGt ^ Q j-QgjQji j\ variatlou in the fragmentation 
scale by a factor of four leads to only a (2-3)% change 
in the total predicted cross section. Similarly a change 
in default set of fragmentation functions (BFG Set 1 to 
BFG Set 2) results in a cross section change of < 1%. 

A possible contribution to the theoretical cross section 
from threshold rcsummation has been estimated (43 | for 
inclusive direct photon production at the Tevatron and 
found to be < (2.5 - 3.0)% for p^ < 350 GcV. 

The experimental systematic uncertainties are reduced 
further by measuring the ratios between the differen- 
tial cross sections D = a / dp^dy'^ dy^'^^ in the differ- 
ent regions. Most of the systematic uncertainties re- 
lated to the identification of central photons then cancel, 
and only systematic uncertainties related to the pp 
7 -l- jet -|- X event purities and the jet selection efficiency 
(when measuring ratios between central and forward jet 
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FIG. 8: Ratio of the predicted cross section with ^.r^fj = 
Prfiv*) to those with Hr^fj = Pt in each measured region. 



regions) remain. Measured ratios between the differen- 
tial cross sections in the different regions are presented in 
Fig. [To] and Tables IVUXl The overall experimental uncer- 
tainty is largest in the first and last bins and ranges 
from (3-9)% across most of the range. The NLO QCD 
predicted cross section ratios estimated using jetphox 
are also presented for scale choices Hr,fj = PTf{y*)i 
lJ-R,Fj = 0.5pJ/(y*), and fiR.F.f = 2pJ/(y*). The scale 
uncertainty of the predicted ratios is < 3% and about 
(3.5-7.5)% for the ratio of cross sections in the two for- 
ward jet rapidity intervals. The shapes of the measured 
ratios between the cross sections in the different regions, 
in general, are qualitatively reproduced by the theory. 
A quantitative difference, however, between theory and 
the measurement is observed for the ratios of the central 



> 



jet regions to the forward 1.5 < \y'°^\ < 2.5, 
region, even after the theoretical scale variation is taken 
into account. The ratio between the two forward jet cross 
sections suggests a scale choice fJ.R,F,f — 2p^f{y*). How- 
ever, the ratios of the central jet regions to the forward 
1.5 < lyj'''! < 2.5, < region suggest a theoretical 

scale closer to ^^r,fj — 0-5p^f{y*). 

In summary, the differential cross section 
d'^a/dp^dy''dy^'^^ for the process pp — > j + jet + X 
is measured for central photons (|y'''| < 1.0) separately 
for four different rapidity configurations between the 
leading photon and the leading jet. The data cover six 
orders of magnitude in the cross section as a function of 
p^ for events with jets in < 0.8, and extend the 

kinematic reach of previous photon plus jet measure- 
ments. Next-to-leading order QCD predictions, using a 
few different modern parameterizations of parton dis- 
tribution functions, are unable to describe the shape of 
the dependence of the cross section across the entire 
measured range. Similarly, theoretical scale variations 
are unable to simultaneously describe the data-to-theory 
ratios in each of the four measured regions. Thus, 



TABLE V: Ratios of 
ly^^'l < 0.8, y-^-y'^' < 
yT.yJ^t > rapidity rej 



the differential cross sections in the 
rapidity region to the \y'°^\ < 0.8, 
;ion. 



7 

Pt 


bin 


Ratio 




<5''syst 




(GeV) 


(r) 


(%) 


(%) 


(%) 


30 - 


- 34 


0.81 


0.4 


9.1 


9.1 


34 - 


- 39 


0.81 


0.4 


7.8 


7.8 


39 - 


- 44 


0.81 


0.6 


6.0 


6.0 


44 - 


- 50 


0.80 


0.7 


4.8 


4.8 


50 - 


- 60 


0.80 


0.8 


3.7 


3.8 


60 - 


- 70 


0.81 


1.3 


3.1 


3.3 


70 - 


- 80 


0.80 


1.8 


2.9 


3.4 


80 - 


- 90 


0.81 


2.5 


2.8 


3.8 


90 - 


110 


0.83 


2.6 


2.8 


3.8 


110 - 


- 130 


0.82 


4.3 


2.7 


5.0 


130 - 


- 150 


0.81 


6.4 


2.4 


6.9 


150 - 


- 170 


0.93 


9.3 


2.2 


9.5 


170 - 


- 200 


0.90 


10.8 


2.1 


11.0 


200 - 


- 230 


0.81 


17.6 


2.0 


17.7 


230 - 


- 300 


0.82 


20.4 


1.8 


20.4 


300 - 


400 


0.67 


58.0 


0.1 


58.0 



TABLE VI: Ratios of the differential cross sections in the 



<0.8, y^-y'°^ > rapidity region to the 
1.5< ly^''*! <2.5, y'^-y'''^ > rapidity region. 



bin 


Ratio 


^J'stat 


Sr syst 




(GeV) 


(r) 


(%) 


(%) 


(%) 


30 - 34 


1.85 


0.4 


11.7 


11.7 


34 - 39 


1.99 


0.5 


10.0 


10.0 


39 - 44 


2.15 


0.7 


8.6 


8.7 


44 - 50 


2.37 


0.9 


7.6 


7.7 


50 - 60 


2.70 


1.0 


6.7 


6.8 


60 - 70 


3.14 


1.6 


5.8 


6.1 


70 - 80 


3.66 


2.5 


5.0 


5.5 


80 - 90 


4.28 


3.6 


4.5 


5.8 


90 - 110 


4.97 


4.0 


4.2 


5.8 


110 - 130 


7.00 


7.5 


3.9 


8.4 


130 - 150 


9.01 


12.6 


3.6 


13.1 



TABLE VII: Ratios of the differential cross sections in the 
|y^°*| <0.8, y'^-y'"*' < rapidity region to the 



1.5<|yj'=*|<2.5, 



y^-y^'^ > rapidity region. 



bin 


Ratio 


Sr Stat 


Sr ayat 


5r:„t 


(GeV) 


(r) 


(%) 


(%) 


(%) 


30 - 34 


1.51 


0.4 


11.9 


11.9 


34 - 39 


1.62 


0.5 


10.6 


10.6 


39 - 44 


1.74 


0.7 


9.1 


9.1 


44 - 50 


1.90 


0.9 


8.0 


8.0 


50 - 60 


2.16 


1.0 


7.0 


7.0 


60 - 70 


2.53 


1.7 


5.9 


6.2 


70 - 80 


2.92 


2.5 


5.0 


5.6 


80 - 90 


3.44 


3.7 


4.5 


5.8 


90 - 110 


4.14 


4.1 


4.3 


5.9 


110 - 130 


5.73 


7.6 


4.0 


8.6 


130 - 150 


7.27 


12.8 


3.7 


13.3 
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FIG. 9; The ratios of the pp ^ y + jet + X cross section with 
just the direct (non-fragmentation) contribution to the total 
(direct+fragmentation) cross section estimated with jetphox 
for each measured region. 



TABLE VIII: Ratios of the differential cross sections in the 



1.5<|2/j'=*|<2.5, 
1.5<|yj'=*|<2.5, 



■y 



jet 
.Jet 



< rapidity region to the 
> rapidity region. 



7 

Pt 


bin 


Ratio 




Sr syst 




(GeV) 


(r) 


(%) 


(%) 


(%) 


30 - 


- 34 


0.49 


0.5 


10.9 


10.9 


34 - 


- 39 


0.50 


0.6 


9.7 


9.7 


39 - 


- 44 


0.49 


0.9 


8.0 


8.0 


44- 


- 50 


0.49 


1.1 


6.7 


6.8 


50 - 


- 60 


0.52 


1.3 


5.6 


5.8 


60 - 


- 70 


0.53 


2.2 


5.0 


5.4 


70 - 


- 80 


0.54 


3.4 


4.8 


5.8 


80 - 


- 90 


0.59 


4.9 


4.8 


6.8 


90 - 


110 


0.54 


5.8 


4.7 


7.4 


110 - 


- 130 


0.64 


10.8 


4.5 


11.7 


130 - 


- 150 


0.63 


18.4 


4.3 


18.9 


150 - 


- 200 


0.97 


32.5 


4.1 


32.8 



the data presented in this Letter, show a need for an 
improved and consistent theoretical description of the 
7+jet production process. 
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TABLE IX: Ratios of the differential cross sections in the 
lyj"*! <0.8, y^-i/^"' > rapidity region to the 
1.5< Ij/-'''*! <2.5, j/'^-yj''* < rapidity region. 



1 
Pt 


bin 


Ratio 




<5''syst 




(GeV) 


(r) 


(%) 


(%) 


(%) 


30 - 


- 34 


3.81 


0.4 


10.4 


10.4 


34 - 


- 39 


4.00 


0.5 


8.8 


8.8 


39 - 


- 44 


4.39 


0.8 


6.9 


6.9 


44 - 


- 50 


4.82 


1.0 


5.5 


5.6 


50 - 


- 60 


5.23 


1.2 


4.6 


4.7 


60 - 


- 70 


5.97 


1.9 


4.3 


4.7 


70 - 


- 80 


6.81 


2.8 


4.5 


5.3 


80 - 


- 90 


7.20 


4.1 


4.6 


6.1 


90 - 


110 


9.21 


4.8 


4.6 


6.7 


110 - 


- 130 


10.91 


8.7 


4.6 


9.9 


130 - 


- 150 


14.31 


14.8 


4.4 


15.4 


150 - 


- 200 


38.29 


23.9 


4.2 


24.2 



TABLE X: Ratios of the differential cross sections in the 



|yj°'| <0.8, y'' ■y''^^ < rapidity region to the 
1.5< <2.5, -y^"^ < rapidity region. 



Pt 


bin 


Ratio 


Sr atat 


Srayat 


5noT 


(GeV) 


(r) 


(%) 


(%) 


(%) 


30 - 


- 34 


3.10 


0.4 


10.7 


10.7 


34 - 


- 39 


3.25 


0.5 


9.5 


9.5 


39 - 


- 44 


3.55 


0.8 


7.4 


7.5 


44 - 


- 50 


3.86 


1.0 


6.0 


6.1 


50 - 


- 60 


4.18 


1.2 


4.9 


5.0 


60 - 


- 70 


4.81 


1.9 


4.5 


4.9 


70 - 


- 80 


5.43 


2.9 


4.5 


5.4 


80 - 


- 90 


5.80 


4.1 


4.6 


6.2 


90 - 


110 


7.67 


4.9 


4.7 


6.8 


110 - 


- 130 


8.93 


8.8 


4.7 


10.0 


130 - 


- 150 


11.55 


14.9 


4.5 


15.5 


150 - 


- 200 


35.48 


23.9 


4.3 


24.3 
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